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ABSTRACT 



Aims. A new method to constrain the cosmological equation of state is proposed by using combined samples of gamma-ray bursts 
(GRBs) and supernovae (SNela). 

Methods. The Chevallier-Polarski-Linder parameterization is adopted for the equation of state in order to find out a realistic approach 
to achieve the deceleration/acceleration transition phase of dark energy models. 

Results. As results, we find that GRBs, calibrated by SNela, could be, at least, good distance indicators capable of discriminating 
cosmological models with respect to ACDM at high redshift. Besides, GRBs+SNela combined redshift-distance diagram puts better 
in evidence the change of slope around redshift z ~ 0.5 which is usually addressed as the "signature" of today observed acceleration. 
This feature could be interpreted, in more standard way, by the red sequence in galaxy clusters. 

Key words. Gamma rays : bursts - Cosmology : cosmological parameters - Cosmology : distance scale 



1. Introduction 

^ From an observational viewpoint, one of the fundamental ques- 
. tion of cosmology is measuring cosmological distances and then 
to build up a suitable and reliable cosmic distance ladder. This is- 
sue has recently become even more important due to the evident 
degeneracy of several dark energy models with ACDM, desp ite 
the advent of the so-called Precision cosmolosv, (lEllis 199 9i). 

More precisely, in the last two decades, a class of accurate 
I standard candles, the Supernovae la (SNela) has been highly 
studied and the results obtained from the use of these objects 
led to the surprising discovery of the appar ent acceleration of 
the cosmic Hubble flow (for a review see ( Kowalski et al 2008^). 
However these objects are hardly detectable at redshifts higher 
' than ~ 1.5, so we need distance indicators at higher redshifts in 
order to remove the disturbing degeneration of dark energy mod- 
els today affecting the current cosmological picture (ACDM is 
a good approximation of the observed Universe, also if there is 
yet no theoretical basis about the nature of its components, but 
the issue of global evolutio n is far from being a ddressed, for 
a comprehensive review see fcopeland et al 2006*)). A possible 
way out for this problem could be found by adopting Gamma 
Ray Bursts (GRBs) as distance indicators also if it is premature, 
at the moment, to speak about standard candles. 

As it is well known, GRBs are the most powerful explosions 
in the Universe: the most likely scenarios for their generation 
are the formation of massive black holes or the coalescence of 
binary stellar systems. These events are observed at considerable 
distances, so there are several efforts to frame them into the stan- 
dard of cosmological distance ladder. In literature, there are sev- 
eral models that give acco unt for the GRB formation. The stan- 
dard model eszaros 2006.) predicts the formation of a black 



hole originated by a massive star whose core is going to collapse. 
Alternatively, the GRB phenomenon could be generated during 
an accretion episode followed by a merging event which gives 
rise to a jet- like outflow. Eor the class of short GRBs (which time 
duration is less than 2 seconds), the candidates a re mergers of 
neutron stars. Another model, (Ruffini et al 2008h . includes dif- 
ferent central energy sources and the formation of charged black 
holes. 

All these scenarios retain essentially a similar shock phe- 
nomenon: a "fireball" or a "fireshell". Here we will not go into 
details, but it is worth stressing that none of these models is in- 
trinsically capable of connecting all the observable quantities. 

Despite of the poor knowledge of the GRB mechanism, 
it seems that GRBs could be used as reliable distance in- 
dicators. In fact there exist several observational correla- 
tions among the photometric and spectral properties of GRBs 
which point out that it c ould be realistic to suppose them 
as distance indicators, (iBasilakos & Perivolaropoulos 20081: 
iGhirlanda et a l 2006). Nevertheless the origin of these spectro- 
scopic and photometrical correlations is not known very well 
and there are several efforts to interpret the behavior of GRB 
feat ures in a coherent \yay, by relatively simple scenarios (e.g. 
see (iDainotti et al 20081: iGhisellini et al2 008)). Succeeding in 
explain the mechanism that generates the GRBs is one of the 
objectives of the modern astrophysics and to clarify these ob- 
served correlations in this context would make GRBs as reliable 
distance indicators. A complete revi ew of the existin g luminosity 
relations for GRBs can be found in (ISchaefer 2007h . 

In this paper, we consider two relations, the one by Liang- 
Zhang ( LZ), (Liang & Zha ng 2005), and the one by Ghirlanda 
(GGL), (iGhirlanda et al 20041) . They are the only 3 -parameters 
relations and have less scatter with respect to the theoretical 
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best fit than the other 2- parameters ones. In a recent paper, 
(ICapozziello & Izzo 20081) . starting from a sample of GRBs, a 
ORB -Hubble diagram has been derived considering these rela- 
tions. However it is worth noticing that the calibration of the 
used relations has been necessary in order to avoid the circular- 
ity problem. This means that all the relations need to be cal- 
ibrated for every set of cosmological parameters. Indeed, all 
ORB distances, obtained in a photometric way, are strictly de- 
pendent on the cosmological parameters since, currently, there 
is no low-redshift (z up to 0.2-0.3) set of GRBs to achieve a 
cosmology-independe nt calibration. In order to overcome this 
difficulty, Liang et al., (Liang et al 2008), proposed a method in 
which several GRB -relations have been calibrated by SNela. In 
fact, supposing that our relations work at every redshift and that, 
at the same redshift, GRBs and SNela have the same luminos- 
ity distance, it becomes possible, in principle, to calibrate the 
GRB -relations using an interpolation algorithm. In this way, it 
becomes possible building a GRB -Hubble diagram by calculat- 
ing the luminosity distance for each GRB with the well-known 
relation between the luminosity distance di and the energy-flux 
ratio of the distance indicators, i.e. 



di- 



(1) 



where Eiso is the isotropic energy emitted in the burst and 5"^^^^ 
is the bolometric fluence corrected to the rest frame of the source 
in consideration. This result can be connected to the Hubble se- 
ries, dVisser 2004). and the density parameters Dm and Da can 
be obtained. The results in ( Capozziello & Izzo 2008) were in 
agreement with t he other observations b ut the estimation of the 
CPL-parameters (IChevallier et al 20011) . which is an advanta- 
geous parameterization of cosmological Equation of State (EoS) 
w = p/p, (see for a brief revie w of the various param eteri- 
zations of the EoS the work by ( Barboza & Alcaniz 2008,) ). is 
only marginally consistent with the data existing in the liter- 
ature, (e.g. see the already cited work by Visser). The reason 
of this disagreement is due to the fact that the method used in 
(ICapozziello & Izzo 20081) works very well at redshift less than 
z ^ 1 while CPL parameters are supposed to work also at high 
redshift. 

The aim of this paper is to take into account a cosmological 
EoS working at any redshift, using GRBs as tracers and adopting 
again the CPL parameterization. The layout of the paper is the 
following: in Sect. 2, we discuss the method which, in principle, 
should allow to obtain a cosmology-independent formulation of 
the luminosity distance and then of the distance modulus. Sect. 3 
is devoted to a discussion of the GRB luminosity-relations con- 
sidered in this work. In Sect. 4, we illustrate the fitting of the 
data obtained by these relations while results and perspectives 
of the approach are discussed in Sect. 5. 

2. The cosmological model 

The goal is to obtain an analytic formulation of the Hubble di- 
agram valid, in principle, at any redshift. Let us start from the 
Friedmann equation 



87rG kc^ 
3 ^ ' 



(2) 



We obtain, by some algebra, the following equation in terms of 
the density parameter 



3(w+l) 



(3) 



where the subscript indicates the present value of the pa- 
rameters. From now onwards, we take into account a spatially 
quasi-flat Universe, ^ 0, the contribution of the curvature 
will be neglig ible and we have Q p ^ 1, as suggested by the 
latest CMBR dKomatsu et al 2008*) and the SNela observations 
dKowalski et al 20081) . However in the final section, we will do 
a test to verify this assumption with observations coming from 
GRBs. Now if we translate in terms of redshift z. 



— = 1 +z, 

a 



the previous equation reduces to 

if2(^)=/f2(l+^)3(w+l) 



(4) 



(5) 



The w-parameter indicates the EoS w = p/p, where p and p 
are the pressure and the matter-energy density of the Universe, 
resp ectively. Consider ing the CPL parameterization of the EoS, 
(Che vallieretal200lh : 



w(z) = Wo + w« - 



and substituting into Eq.(|5]), we obtain: 



H(z) = Ho 



(l+z)^^-^-^^^exp(-^ 
\2(1 +z) 



(6) 



(7) 



which enters directly in the expression of the distance modulus 
Kz) = -5 + 5 log d,(z), (8) 
where di(z) = c(l + z)Di(z) and where 



D,{z) 
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This means that an analytic expression for yu can be achieved. 
The integral Di in Eq.(|9]) can be solved giving a Gamma function 
of the first kindQ: 



exp I IF 



1 -h 3wo + 3w, 



3Wa 



2(1+ a 



1 +31^0+31^0 
2 

^=0 



(10) 



Substituting such an expression in the distance modulus, we 
obtain a model for data fitting which could work, in principle, at 
any z. It is important to stress that the obtained expression for the 
Hubble parameter H(z) is independent of the density parameters, 
Qm and Q\, provided that their sum is equal to L 

It is worth noticing that we are using the CPL parameteri- 
zation not only for the dark energy component, but for the total 
energy-matter density of the Universe. This assumption works 
because dark and baryonic matter are contributing with a null 
pressure while the radiation component is negligible in matter- 
and dark energy-dominated eras. Furthermore, the analytical for- 
mulation which we are adopting for the luminosity distance is 
assumed valid at any redshift z. 



^ In our case, the variable of the Gamma function, z, is always posi- 
tive so that we have no problem of discontinuity in applying the Gamma 
function in the following calculations. 
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3. GRBs luminosity relations 

As we said, GRBs are the most powerful explosions in the 
Universe, so they can be observed up to very far distances. 
Theoretically, they can be observed up to a redshift of the order 
z ~ 10. Some considerations on the spectrum of such sources 
are due at this point. The electromagnetic emission spectrum 
of GRBs ranges from radio up to gamma wavelengths, but the 
main bulk of emission is in the gamma band. In the last years, 
thanks to several spacecraft missions capable of observing this 
high energy region, the main features of GRBs have been bet- 
ter known. Recently, some photometric and spectroscopic rela- 
tions between GRB observables have been found and then the 
hypothesis that these objects could be considered suitable dis- 
tance indicators has been seriously considered. Nevertheless, up 
to now, there is no theoretical model that fully explains these re- 
lations so the GRBs cannot be considered as standard candles in 
a proper sense. For a de tailed review of the observational fea- 
tures see (ISchaefer 2QQ7h . 

Here, we are taking into account the existing 3 -parameter 
relations. This choice has been done because these relations 
put the better constraints on the data giving less scatter be- 
twe en the the oretical relation and the experimental data (e.g. 
see (ISchaefer 2007D). The first relatio n is the so-called Liang- 
Zhang relation. (iLia ng & Zhang 2QQ5I) . which allows to connect 
the GRB peak energy, Ep, with the isotropic energy released in 
the burst, Eiso, and with the jet break - time of the afterglow op- 
tical light curve in the rest frame, measured in days, t^, that is 



XogEiso =a + bi\og 



?>mkeV 



+ ^2 log 



tb 



(1 -h z)lday 



(11) 



where a and Z?/, with / = 1,2, are calibration constants. 

The other one is the relation given by Ghirlanda et al. 
(iGhirlanda et al 20 d3). It connects the peak energy Ep with the 
collimation-corrected energy, or the energy release of a GRB jet. 



Ey, where 



(12) 



and Fheam = 1 - COS 0, with Ojet the jet opening angle defined in 
(ISari et al. 1999h : 



Ojet = 0.163 



3/8/ noTJy 
Zl \ Eiso,52 
^52 



1/8 



(13) 



where Eiso^si = Eisoll^ ergs, is the circumburst particle 
density in 1 cm"^, and riy the radiative efficiency. The Ghirlanda 
et al. relation is 



log Ey = a-\- b log 



300y^^y 



(14) 



where a and b are two calibration constants. 

From these relations, we can obtain directly the luminosity 
distance di from the well-known formula which connects di with 
the isotropic energy Eiso and the bolometric fluence S boh • 



di = 



47iS 



bolo 



(15) 



from which it is easy to compute, for each GRB, th e distance 
modulus yu = and its error given by dLiang et al 2008|) : 



{2.5criogEi,of + '^^^CTS tola bolo f 



(16) 



Table 1. Results of the fits. SNela is only for the Supernova la 
data, LZ is for the GRBs data obtained from the Liang-Zhang 
relation, GGL for the Ghirlanda et al. one. Note the improvement 
on the w-parameter using the GRBs data in addition to the SNela 
data corrected for the 3 "outlier" GRBs. 



Relation 


Wo 




R' 


SNela 


-0.910 + 0.070 


0.755 + 0.054 


0.983 


LZ 


-1.39 + 0.38 


1.18 + 0.37 


0.817 


GGL 


-1.46 + 0.38 


1.36 + 0.32 


0.812 


LZ + SNela 


-1.15 + 0.10 


0.93 + 0.11 


0.933 


GGL + SNela 


-1.42 + 0.12 


1.24 + 0.13 


0.920 



with criogEi.o ^Sboio obtained from the error propagation ap- 
plied to Eq.JTTI) and Eq.([T4l). Moreover, we assume that the error 
in the determination of the redshift z is negligible, as well as for 
the radiative efficiency r]y. We note also that the assumption of 
a well-known no is a strong hypothesis since the goodness of 
the fits depends, in particular, on this parameter. However, we 
are going to consider the no values for each GRB given in the 
Table[5]also if it lacks a complete and clear phys ical basis for the 
considered relations, (iFriedman & Bloom 20051) . The GRB data 
sample is taken from the already cited work by Schaefer. We 
take into account 27 events wi th extremely precise data. Such a 
sample is the same adopted in dCapozziello & Izzo 2008h . 

4. The data fitting 

The next step is the fit of the GRB sample with the empirical re- 
lations, Eqs.([TT1),([T4l). described in Sect. 3. The aim is to achieve 
an estimate of the GPL parameters and consequently to deter- 
mine the trend of the EoS at any redshift, using the analytical 
relation, Eq.JTOt. As we said, we are considering the same sam- 
ple of 27 GRBs used in ( Capozziello & Izzo 2008), see Table 
[51 in which we hav e added the sample of SNela by the Union 
Supernova Survey, (iKowalski et al 20081) . 

The numerical results of the fits are shown in Table [T] where 
we obtain a robust estimation of the CPL parameters for both 
the relations used, with and without SNela data. An immediate 
comparison is done with the best fit applied only to the SNela 
sample. It is evident how adding GRBs to SNela data improves 
the knowledge and the precision on the EoS parameter w. In Fig. 
[U the best fit curve, in the case of LZ relation, is plotted while, 
in Fig. [21 it is plotted the best fit curve for the SNela data using 
the theoretical model described previously. 



In order to measure the goodness of the fit, we use the test 
for an accurate reliability, see Table [H The R^ test is a measure 
of how successf ul the fit is in explainin g the variation of the data 
(see for details ([Draper & Smith 1998[) ). Exactly an R^ close to 
1.0 indicates that we have accounted for almost all of the vari- 
ability with the data specified in the model. As standard, the R^ 
test is the square of the correlation between the response values 
and the predicted response values, that is: 



R' 



1 



SSE 
SST 



1 



i:Um(yi-yi)' 



(17) 



where SSE is the sum of the squares due to errors and it mea- 
sures the total deviation of the response values from the fit and 
SST is the sum of squares about the mean: y is the predicted 
response value, y is the mean value and the w/ are the weights on 
the values. 
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Fig. 1. Redshift-Distance modulus diagram for the GRB+SNela 
sample. The black dots are the GRBs , the blue ones are the 
SNela. The red line is the best fit obtained from the data, with the 
dashed line representing the confidence limits at 3cr. The error 
bars on the Supernova data are not represented because they are 
negligible. 
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Fig. 2. Redshift-Distance modulus diagram for the SNela sam- 
ple only. The black dots are the Supernova data while the red 
line is the best fit obtained from the data. The error bars on the 
Supernova data are not represented because they are negligible. 



A further consistency test for the adopted samples of GRBs 
and SNela can be derived considering the relation (log J/)^/"^ - 
red shift vs z. In such a way, data from GRBs and SNela can 
be better separated in order to track the whole trend of both sets 
and showing, in particular, the spreading of GRB distribution. In 
Fig. [3l the results of this relation is plotted and it is reasonable 
to conclude that SNela could be used to calibrate GRBs at lower 
redshifts. In particular, it is evident, even considering the only 
SNela data, a change in the trend at redshift included between 
z = 0.3andz = 0.5. 

In particular, the extension of the Supernova Hubble 
Diagram with the GRB data can be used to improve our knowl- 
edge of the trend at high redshift. In this way, using also the GRB 
data, we build a plot. Fig. [5j where the distance modulus ji ver- 
sus the redshift z, in a logarithmic scale, is plotted. The best fit 
curve, obtained with the method described previously, Eq.dTOl), 
is also reported. However in order to better point out this varia- 




Fig. 3. Comparison between GRBs- SNela data using the relation 
(logJ/)^/^ for each data sample. Blue dots are the SNela while 
the red star are the GRBs. We are considering redshifts between 
and 3, in order to better analyze the data trend in the overlap- 
ping redshift range. 



Fit Residuals 




Fig. 4. Comparison between the best fit of fi and the observed 
distance modulus jiobs at any redshift. The black dots are the 
GRBs data and the red line is the best fit curve representing the 
theoretical distance modulus. 

tion, we have done an extended analysis. Fig. [6l where we con- 
sidered the GRB-^Sn sample up to a certain value of the redshift 
Zt and its complementary, and then we fitted these samples with a 
simple curve of the type a/ -h/?/ log z, with the cut-redshift Zt rang- 
ing from Zt = OA to Zt = 0.6 with redshift-step 0.1. The result of 
this analysis is given in Table [21 where we note that at redshift 
greater than Zt = 0.5, the two best fit curves tend slowly to coin- 
cide, suggesting that something happens at redshift smaller than 
this value. This result led in the past to the conclusion that the 
universe is accelerating. 



It is suggestive to link the observed behaviour of the comso- 
logical accelletration to the the formation and evolution of astro- 
phy sical structures. F or instance, the infalling rate of field galax- 
ies dStott et al. 2009b could be related to the overal expansion 
rate. 
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Fig. 6. Redshift-Distance modulus diagram for the GRB+SNela samples versus log z as described in the text. Note the "linearization" 
of the two best fit curves of the two data samples, done with the model at + hi log z, obtained considering only ORB and Sn data up 
to and beyond, respectively, the cut-redshift value Zt, with Zt ranging from Zmin =0.1 and Zmax = 0-6 with step of Zstep = 0.1. 

Table 2. Results of the logarithmic fit described in the text and plotted in fig. [6] 



Zt 


ai 


bi 


a2 


bi 


0.1 


44.07 


2.52 


43.48 


2.25 


0.2 


44.07 


2.52 


43.55 


2.28 


0.3 


44.08 


2.55 


43.81 


2.34 


0.4 


44.07 


2.65 


43.85 


2.35 


0.5 


44.06 


2.48 


43.90 


2.37 


0.6 


44.06 


2.52 


43.95 


2.38 



It would be very interesting to asses quantitatively such a 
relation, from the massive galaxy clusters with some good in- 
dicator at that redshift in order to clarify a possible correlation 
between these two phenomena. All of these cosmological and 
astrophysical hints lead us to reconsider from a diff'erent point 
of view the possible cause of the apparent acceleration of the 
Universe which could be addressed to some more astrophysical 
standard eff'ects, (Izzo et al. in preparation). 

In Fig. m it is plotted the comparison between the theoret- 
ical j^th and the observed distance modulus j^obs at any redshift, 
the residual plot. A smooth trend up to z ~ 3.5 in the resid- 
ual curve can be immediately detected. Beyond this limit, we 
have 3 GRBs that exceed, by the same side, the 3cr confidence 
limit of the best fit. This discrepancy is clear in Fig. [5j where 
is plotted the best fit for the combined sample in the case of 
LZ relation with a logarithmic scale for the redshift. It is worth 
noticing that a similar, but opposite discrepancy was obtained by 



(I Perivolaropoulos & Shafieloo 2008h using only the SNela sam- 
ple. In that paper, the authors found a brightening for the SNela 
at high redshift, while here, we find a sort of "darkening" for the 
GRBs at redshift higher than z = 4. 

This fact is fundamental for the goodness of the fit because 
these GRBs represent the most distant objects that one can use 
to make such an analysis and their weight on the fit is very high, 
in the sense that they appear to be not accurate distance indica- 
tors. There could be several explanation for this anomalous GRB 
brightness at high redshift and the most likely are the following: 



- there is some process of absorption of gamma radiation 
where the GRB y photons may interact with the very low en- 
ergy photons in coming from the cosmic t hermal background 
radiation, (Zd ziarski and Svensson 1989b : 

- these 3 GRBs could be outliers; 
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Redshiftz 

Fig. 5. Redshift-Distance modulus diagram for the GRB+SNela 
sample versus redshift in logarithmic scale. 

- the Circum Burst Density (CBD) in the host galaxy at that 
epoch, z > 4, could be different from the CBD in galaxies at 
low redshift; 

- there could be some high-energy bias: since at that dis- 
tances only very powerful GRBs can be observed, some 
high-energy process, involving very energetic y-photons, 
(Kelner et al. 2008^ iRazzaque et al. 200 9), could happen so 
that the flux received by our detectors is dimmed; 

- it could be that the CPL parameterization, or the ACDM 
model, is a bad approximation for the cosmological EoS. 

It is very interesting to note that this phenomenon is very 
similar to the change in the trend of the Supernovae la data, as 
we mentioned earlier. However, at these redshifts is very difficult 
to make any kind of astrophysical constraint, since there is not 
enough information about this region of the Universe. 

Nevertheless, we repeat the analysis described above with- 
out these 3 GRBs obtaining a better value than the previous one 
for the test. The results of these corrected fits are shown in 
the table m In Eig. [71 it is plotted the best fit with this corrected 
sample. Erom these results we conclude that the complete sam- 
ple gives different results from the corrected sample, the first one 
suggesting a phantom/quintessence regime for the present epoch 
while the second one can be enclosed in the case of an accelerat- 
ing ACDM model. This last result is confirmed by the following 
analysis, where we have performed a Monte-Carlo-like proce- 
dure for the comparison of the results with the usual likelihood 
estimator given by 



X 



(j^thiZi) - JdobsiZi))^ 



(18) 



in the context of a ACDM model of the Universe, where jith is 
the distance modulus computed from the Eq.© and Eq.(|9]), Zi is 
the observed redshift for each GRB and ct/ the observed distance 
modulus uncertainty. The results of this analysis are shown in the 
TableO where we can see the improvement obtained by the GRB 
sample corrected for the 3 "wrong" GRBs. In Eig.®, it is shown 
the contour plot of the corrected sample, where the boundaries 
correspond to Icr, 2cr and 3cr confidence levels. In this, case we 
do not consider a purely fiat geometry, as we previously said in 
the Sect. 2, but we constraint the /:-parameter of Eq.© to vary 
between the value -0.05 < k < 0.05 so that we take account of 
a possible small contribution due to the curvature density. 




1.5 2 
Redshift z 



Fig. 7. Redshift-Distance modulus diagram for the corrected 
GRB-hSNela sample. The red line is the best fit obtained from 
the data, with the dashed line representing the confidence bound- 
ing at 3(T. 



We have adopted a similar procedure in the case of an EoS 
evolving with redshift and where jith is obtained by the Eq.ifTOb. 
The result of this analysis is plotted in Eig. [9] where the best 
fit value, the cross in the figure, corresponds to the value wq = 
-0.84 ± 0.14 and Wa = 0.72 ± 0.06, in a good agreement with 
the results obtained, see Table IH using our theoretical relation, 
Eq.dTOl). 

Summarizing, from this analysis, we conclude that the cor- 
rected sample agrees fairly well with the ACDM model with 
a mall contribution of the curvature parameter, being equal to 
k = 0.01 ± 0.04. In other words, the method delineated in the 
Sect.2 seems a good approximation of the observed cosmogra- 
phy and agrees very well with the ACDM model so that we can 
argue that GRBs could be good distance indicators at redshift 
values up to z = 4. 




Fig. 8. 68%, 95% and 98% constraints on and Qa, see Eig.( 
[8]) obtained from UNION sample and the GRB sample corrected 
for the 3 wrong GRBs. 
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Fig. 9. 68%, 95% and 98% constraints on wq and Wa obtained 
from UNION sample and the GRB sample corrected for the 3 
wrong GRBs. The cross represents the best fit value and it is in 
a good agreement with what found using the theoretical model 
described in Sect. 2. 



5. Discussion and Conclusions 

Starting from the Friedmann equation, we have investigated a 
new method to constrain the cosmological Equation of State at 
high redshifts. The working hypothesis lies on the use of GRBs 
as distance indicators at high redshift, well beyond the distance 
where SNela are detected up to date. The CPL parameteriza- 
tion for the EoS has been explicitly used for the whole matter- 
energy content of the Universe as a suitable approach to investi- 
gate the parameter w = w(z) and discriminate with respect to the 
ACDM model. In particular, regarding the Friedmann equations, 
for which 



(1 + ^) 



(19) 



where q is the deceleration parameter and where the prime 
denotes the derivative with respect to the redshift, we have ob- 
tained, in the case of the LZ relation, with a reliable confidence 
level, the epoch for the transition between the deceleration- 
acceleration phases at a redshift value of z ~ 5. This is a value 
that, also if higher than the redshift of the farther GRB used, 
could be in agreement with current quasar formation scenar- 
ios. Besides, we are in good agreement with the observed phan- 
tom/quintessence regime at present epoch, that is for z ^ 0, we 
obtain w < - 1 . However we have found an anomaly at z ~ 4 and 
beyond for which we rejected 3 GRBs at that distances. Several 
explanations are possible for this problem which can be summa- 
rized as: 

- the eff'ective cosmology is not described by the ACDM 
model 

- the GRBs evolve with the redshift 

- there is some nuclear or electromagnetic processes between 
the y photons and the baryons involved in a GRB that would 
dim the observed flux or fluence. 



Table 4. Results of the fits corrected for the 3 "wrong" GRBs. 
SNela is just for the Supernova la data, LZ is for the GRBs data 
obtained from the Liang-Zhang relation, GGL for the Ghirlanda 
et al. one. 



Relation 


Wq 






LZ + SNela 


-0.95 + 0.01 


0.74 + 0.01 


0.999 


GGL + SNela 


-0.865 + 0.005 


0.66 + 0.005 


0.999 



Nevertheless such issues still remain unsolved and it repre- 
sents a challenge. With this fact in mind, we have performed 
the same analysis without these 3 GRBs obtaining diff'erent re- 
sults from the previous ones. In particular, we rejected the today 
phantom regime by this new analysis, obtaining for wq a value in 
agreement with the ACDM model at present epoch. The method, 
also if preliminary, seems to indicate that GRBs could be actu- 
ally used as standard candles once a reliable unified model for 
the photometric and spectroscopic quantities is achieved (reli- 
able results in this sense are presented in ( Ghisellini et al 2008) ). 
However, more robust samples of data are needed and more real- 
istic EoS (with respect to the simple perfect fluid models) should 
be taken into account i n order to suitably trac k redshift at any 
epoch (see for example (ICapozziello et al. 2006)). 

With the improving of the observations, in particular with the 
launch of new satellites devoted to the GRB surveys, as Fermi- 
GLASlElandAGIL^I, one should be able to expand the samples 
of GRBs, possibly with data coming from objects at higher red- 
shift. 

In summary, considering these preliminary results, it seems 
that GRBs could be considered as a useful tool to remove de- 
generation and constrain self-consistent cosmological models. 
Furthermore the matching with other distance indicators would 
improve the consistency of the Hubble distance-redshift diagram 
by extending it up to redshift 6-7 and over. 

We wish to thank Riccardo Benini for useful discussions and 
suggestions. 
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Table 3. Cosmological density parameters, with uncertainties computed at Icr confidence Hmit, obtained from the MC-Hke proce- 
dure 



Sample 






Ha 






UNION H 


hGRB 


0.26 + 0.14 


0.73 + 0.14 


0.01 + 0.04 


1.032 


UNION H 


h ORB corrected 


0.25 + 0.10 


0.74 + 0.135 


0.01 + 0.035 


1.00027 



Table 5. GRBs Data Sample 



GRB 


z 


Ep (keV) 


Sboio (erg cm ^) 


tjet (days) 


Ojet (deg.) 


no {cm ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


970508 


0.84 


389 + 40 


8.09E-6 + 8.1E-7 


25 + 5 


23 + 3 


3.0 + 2.4 


970828 


0.96 


298 + 30 


1.23E-4 + 1.2E-5 


2.2 + 0.4 


5.91 + 0.79 


3.0 + 2.4 


980703 


0.97 


254 + 25 


2.83E-5 + 2.9E-6 


3.4 + 0.5 


11.02 + 0.8 


28.0 + 10 


990123 


1.61 


604 + 60 


3.11E-4 + 3.1E-5 


2.04 + 0.46 


3.98 + 0.57 


3.0 + 2.4 


990510 


1.62 


126 + 10 


2.85E-5 + 2.9E-6 


1.6 + 0.2 


3.74 + 0.28 


0.29 + 0.14 


990705 


0.84 


189 + 15 


1.34E-4 + 1.5E-5 


1 + 0.2 


4.78 + 0.66 


3.0 + 2.4 


990712 


0.43 


65 + 10 


1.19E-5 + 6.2E-7 


1.6 + 0.2 


9.47 + 1.2 


3.0 + 2.4 


991216 


1.02 


318 + 30 


2.48E-4 + 2.5E-5 


1.2 + 0.4 


4.44 + 0.7 


4.7 + 2.8 


010222 


1.48 


309 + 12 


2.45E-4 + 9.1E-6 


0.93 + 0.1 


3.03 + 0.14 


3.0 + 2.4 


011211 


2.14 


59 + 8 


9.20E-6 + 9.5E-7 


1.56 + 0.16 


5.38 + 0.66 


3.0 + 2.4 


020124 


3.20 


87 + 18 


1.14E-5 + l.lE-6 


3 + 0.4 


5.07 + 0.64 


3.0 + 2.4 


020405 


0.70 


364 + 90 


l.lOE-4 + 2.1E-6 


1.67 + 0.52 


6.27 + 1.03 


3.0 + 2.4 


020813 


1.25 


142 + 14 


1.59E-4 + 2.9E-6 


0.43 + 0.06 


2.8 + 0.36 


3.0 + 2.4 


021004 


2.32 


80 + 53 


3.61E-6 + 8.6E-7 


4.74 + 0.5 


8.47 + 1.06 


30.0 + 27.0 


030226 


1.98 


97 + 27 


8.33E-6 + 9.8E-7 


1.04 + 0.12 


4.71 + 0.58 


3.0 + 2.4 


030328 


1.52 


126 + 14 


6.14E-5 + 2.4E-6 


0.8 + 0.1 


3.58 + 0.45 


3.0 + 2.4 


030329 


0.17 


67.9 + 2.3 


2.31E-4 + 2.0E-6 


0.5 + 0.1 


5.69 + 0.5 


1.0 + 0.11 


030429 


2.66 


35 + 12 


1.13E-6+ 1.9E-7 


1.77 + 1.0 


6.3 + 1.52 


3.0 + 2.4 


041006 


0.71 


63 + 12 


1.75E-5 + 1.8E-6 


0.16 + 0.04 


2.79 + 0.41 


3.0 + 2.4 


050318 


1.44 


47 + 15 


3.46E-6 + 3.5E-7 


0.21 + 0.07 


3.65 + 0.5 


3.0 + 2.4 


050505 


4.27 


70 + 23 


6.20E-6 + 8.5E-7 


0.21 + 0.04 


3.0 + 0.8 


3.0 + 2.4 


050525 


0.61 


81.2 + 1.4 


2.59E-5 + 1.3E-6 


0.28 + 0.12 


4.04 + 0.8 


3.0 + 2.4 


050904 


6.29 


436 + 200 


2.0E-5 + 2E-6 


2.6 + 1 


8 + 1 


3.0 + 2.4 


051022 


0.80 


510 + 22 


3.40E-4+ 1.2E-5 


2.9 + 0.2 


4.4 + 0.1 


3.0 + 2.4 


060124 


2.30 


237 + 76 


3.37E-5 + 3.4E-6 


1.2 + 


3.72 + 0.15 


3.0 + 2.4 


060210 


3.91 


149 + 35 


1.94E-5 + 1.2E-6 


0.33 + 0.08 


1.9 + 0.17 


3.0 + 2.4 


060526 


3.21 


25 + 5 


1.17E-6+ 1.7E-7 


1.27 + 0.35 


4.7 + 1 


3.0 + 2.4 
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